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“Slow” methods = Controlled or 
“dialable” accuracy 


e What do we mean by this? 
— Specify number of digits of accuracy in advance 


— Computer code reliably meets that specification 
(or tells user that 1t cannot) 


* Byproduct: estimate of accuracy of results 
— Error bars 


| 
| Georgia ins “itiutc 


| — | 


Controlled accuracy 


* Few applications demand "high" accuracy 
— Limited resolution of actual objects 
— Material parameters not well characterized 
* Ability to obtain relatively high accuracy 
necessary for controlled accuracy 
— Controlled accuracy does not imply high accuracy 
e Ultimately, controlled accuracy techniques should 
be the most efficient 
— They may be “slow” for now, but... 
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What are the necessary ingredients? 


e Robust formulation 

e High resolution model of object geometry/materials 
* High order representations of currents/fields 

e Treatment of edge/corner/tip singularities 

* Accurate Green's function computations (IE) 

* Means of estimating error 


e Adaptive refinement (h-, p-, hp-) for efficiency 


| 


Typical practice today 


echniques based on “low order” representations 


— Example: RWG basis functions, flat faceted model 
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Whats wrong with RWG? 
Low-order accuracy 


e RWG currents behave as O(h) or O(h?) in the best case 


e RWG Galerkin EFIE SCS behaves as O(h°), a 
superconvergent rate, in the best case 


* [ow order methods make it difficult to improve accuracy 
very much, also difficult to estimate accuracy 
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Opposite of low order: High order 


* Basis functions of polynomial degree p 


— Complete/mixed-order/interpolatory/hierarchical 


* Geometry represented by curved cells 
— Exact/approximate nuum 
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* Computational precision 
— Must be adapted as needed ; 
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High order bases offer 2 advantages 


e Better accuracy (but not necessarily 'dialable") 


e Reduced density of unknowns 
— RWG: 200 unknowns per sq. wavelength 
— WIPL: 30 unknowns per sq. wavelength 
— Recent paper by Ludick, Tonder, Jakobus (EMSS), 


ICEAA 2013: 
Method: RWG HOBF 
Time for 
matrix setup 1.9 hr | 38.6 min 


‘Time for 
matrix solution 7.3 hr 4.5 min 
é 1 


Memory usage 
umo (whole solution) | 26.2 GB | 1.2 GB 
OU nology 








How do we measure “error”? 


e Error estimator used for research purposes: 
— Solve an over-determined system (2N by N) 
— Residual error available as byproduct: 
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— Good correlation with actual error 1n problems where 
we can compute the actual error 
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How do we measure “error”? 


— Residual error defined locally or globally 
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— Slope of residual error curves often useful 
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Slope 
“4 logs (N,) - log, QV) 


Simple (simplest?) example 


* Circular PEC cylinder, radius a = | wavelength 
— IM-to-z MFIE 
— Exact circular model 
— Uniform plane wave excitation 


— MoM solution w/ q Legendre polynomial basis 
functions per cell and point testing at 2q locations 
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Example: Circular cylinder 


MEIE-TM on a 2-D circular cylinder; a=1.0 
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Normalized Residual Error - L2 
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Observation 


e When we see convergence rates (slopes) that are 
the same as the circular cylinder problem, have 
high confidence in the results 

— things are working OK! 


* Reality: don't see this kind of convergence 
without doing (almost) everything right 


What we learned: smooth targets 
(spheres, toroids, prolate spheroid) 


* High order techniques work well for smooth 
targets 
— Slopes of error curves increase with p 
— Observe consistent error slopes 
— More efficient use of computation 
e Biggest challenge: maintaining high precision in 
Green's function calculations for integral eqns 
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What we learned: linear dipole 


* Need accurate Green's function integrations 
* Cannot reduce residual error uniformly without 
including charge singularities at discontinuities 
— special basis for singularity needed for end of hollow 
tube, flat end cap models 
* Helps to have variable cell sizes 


— smaller cells where more rapid variation in excitation or 
current density encountered 
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Current/charge modeling at edges 


* Polynomial basis functions do not properly model 
the current and charge density near the edge 


* Do not obtain high order behavior (steeper slopes) 


e Local residual error does not decrease near surface 
discontinuity 
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High order approach for edges: 
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Consider multiple terms from the wedge solution in 
order to produce residuals with steeper slopes 


Example: 60 degree corner 


Table I. 
Example: Exponents of the degrees of freedom used for a given order, for cells 
near a 60 degree corner (TM case). The corner cells involve twice as many 
degrees of freedom as the other cells, and are twice as large. 


Order of Exponents used in non- | Exponents used in corner cells 
representation q corner cells 


0 -2/5, 0 
-2/5, 0, 1/5, 1 
-2/5, 0, 1/5, 1, 4/5, 2 


0,1,2, 3, 4 —2/5, 0, 1/5, 1, 4/5, 2, 7/5, 3, 
8/5, 4 





-2/5, 0, 1/5, 1, 4/5, 2, 7/5, 3 


Triangular cylinder target 


e TM PEC cylinder, equilateral triangle cross section 
— 7 wavelengths per side 
— Corner-on incident uniform plane wave E!“ کے‎ 
— MFIE "ws < 


e Expansions of order g = 2, 3,4, with g unknowns 
(degrees of freedom) per non-corner cell 


* Add q singular terms in corner cells 


* Uniform cell sizes except corner cells twice as large 
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Polynomial vs. Singular expansions 


q-2 
q = 3 


q = 4 
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Square Cylinder 


— 7 wavelengths per side 
— Corner-on incident uniform plane wave 
— MFIE solution 


TM PEC cylinder, square cross section اع‎ S 


Employ expansions of order g = 2, ...,6, with q 
unknowns (degrees of freedom) per non-corner cell 


Add q singular terms in corner cells 


Uniform cell sizes except corner cells twice as large 
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Observation 


* Convergence rates (slopes) for triangle and square 
cylinder are the same as the circular cylinder 
— true high order behavior 


— things are working OK! 


e Current density for q = 7 has converged to at least 
5 decimal places 


مو مج وج جج ٭ ٭ 9و و9 و ٭ و > و 9> و 
structure 6$90909099?090909990922909296‏ 43-1 
LCN results for TE excitation on one strip‏ 

of an infinite 3-fin of 7 wavelengths width/strip. 
The junction is at cell 1 and the edge is at cell 26. 
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中 让 让 


e local NRE for: 
— poly only 
— poly + end 


poly + junction I 0909000090009000009000000920 
aaAasAaAAaAAaAAAAAAAAAAAAAASA 


poly + end + junction 


wocooooooooocoooococoo2o9o9co93。 
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Recommendation (2D): 


e Starting point for analysis of problems with edges/corners: 


— Use uniform cell sizes, except that corner cells are twice the size of 
other cells 


— Use uniform order g in non-corner cells 
— Use an additional g singular terms in corner cells 


e h-refinement should yield high order behavior 


e Ultimately, employ some form of adaptive h/p refinement 
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Tip/Corner Singularities in 3D 


* Vertex of plane angular sector has exponents associated 
with Lame functions 


— corners of plate and aperture 
— (Satterwhite & Kouyoumjian, 1970) 


* Tip of circular cone has exponents associated with the 
zeros of Legendre functions 


— (Van Bladel, various) 


e Tip of elliptic cone has exponents associated with Lame 
functions 


— (Boersma & Jansen, 1990) 
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3D target 71: rectangular plate 


e Think we know how to handle plate edges based on the 
approach that works in 2D 


— Adapt preceding approach to a surface, determine how much larger 
to make edge cells to balance local residual errors, etc. 


* Need to investigate plate corners 


— Can we reduce residual errors without modeling the corner 
singularity? 
。 Incorporate corner behavior into special basis functions 


— Goal is to produce residual errors that exhibit same slopes as 
circular cylinder problem 
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Rectangular plate corner current 
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Conclusions 


e Today, a wide range of 2D problems can be solved to high 
accuracy even when edge singularities are present 


— have yet to demonstrate "controlled" accuracy with integral 
equations, but we are close 


— need reliable & efficient error estimators, adaptive refinement 
schemes 


e Still have quite a way to go to achieve “controlled” 
accuracy for general 3D problems 


— missing pieces: tip singularities, error estimators, better grasp of 
model error impact, ... 
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